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SATELLITE INFRARED SOUNDINGS FROM NOAA SPACECRAFT 


L.°M. McMillin, D..O. Wark, J. M..Siomkajlo, 
PP. G. Abel, A. Werbowetzki, L. A. Lauritson, 
J. A. Pritchard, D. S. Crosby, H. M. Woolf, 
R. C. Luebbe, M. P. Weinreb, H. E. Fleming, 
F., E. Bittner and C. M. Hayden 
National Environmental Satellite Service 
National Oceanic and Atmospheric Administration 
Washington, D.C. 


ABSTRACT. Data are currently being received from 

a Vertical Temperature Profile Radiometer (VTPR) 
aboard the NOAA 2 spacecraft to produce operational 
atmospheric soundings of temperature and humidity 
on a global scale. This report describes the 

VTPR instrument, its calibration, the procedure 

to obtain "clear radiances" from cloud-contaminated 
radiance measurements, retrieval techniques used 

to obtain temperature and humidity profiles from 
"Clear radiances", the quality checks performed 

on these profiles, and the various forms in which 
data are available to potential users. 


i... INTRODUCTION 


Beginning with the NOAA 2 satellite, launched 15 October 1972, 
satellites of this series will carry instruments to make routine 
observations from which atmospheric temperature soundings are 
derived. This is an operational system designed to provide glo- 
bal upper air data to weather services. The purpose of this 
report is to describe the means by which the operational product 
is derived from satellite measurements. Subsequent chapters 
describe the instruments, ground support equipment, data systems, 
mathematical developments, and the forms of final products. 


Basic concepts of indirect soundings from satellites have been 
tested on the Nimbus 3 and 4 satellites. A Satellite Infra-Red 
Spectrometer (SIRS) (Wark 1970, Wark and Hilleary 1969, and Wark 
et al. 1970) and an Infra-Red Interferometer Spectrometer (IRIS) 
(Conrath et al. 1970, Hanel and Conrath 1969, and Hanel et al. 
1972) were carried on each of these satellites. Output signals 
from these instruments provided measurements of the thermal radi- 
ation emitted by the earth's atmosphere and surface. These 
measurements can be related to the vertical structure of 
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temperature and humidity. A full description of the principles of 
measurement and data reduction is given by Fritz et al. (1972). 


Following initial verification of the reliability with which 
temperature profiles can be reproduced, the results from SIRS 
data were provided on a timely basis to the U.S. National Weather 
Service and were sent through normal telecommunications to other 
weather services (Smith et al. 1972). This quasi-operational 
procedure waS commenced in May 1969 and continued with 
occasional interruptions until the implementation of the fully 
operational product from NOAA 2. 


In anticipation of the increasing demand for quantitative wea- 
ther observations from satellites, the NASA and the NOAA have 
planned and produced a system to obtain operational soundings 
covering most areas of the globe twice daily. Plans have specif- 
ically dealt with the main deterrent to accurate soundings - 
clouds. As will be shown, the effect of clouds can be overcome 
by conducting many measurements within the area from which a 
single sounding is derived. Soundings are obtained from a num- 
ber of sets of measurements; SIRS instruments obtained only 
Single sets of measurements from which soundings were derived. 

A statistical technique (Smith et al. 1972) permits one to deduce 
the emitted radiation from cloud-free areas. 


Each satellite in the NOAA series beginning with NOAA 2 will 
carry Vertical Temperature Profile Radiometers (VTPR's). : 
Duplication of instruments on each spacecraft will assure con- 
tinued operation in the event one fails. The spacecraft are 
described by Schwalb (1972). 


NOAA satellites orbit the earth in 78.3° retrograde orbits each 
115 minutes at an altitude of 1464 km. Figure 1 shows the earth 
projection of seven orbits during a 13-hour period. Solid lines 
indicate the north to south, or "descending", portions of the 
orbits, which occur over mostly sunlit areas; dashed lines are 
the nighttime "ascending" portions of the same orbits. Equator 
crossings of NOAA 2 occur at 9 a.m. (descending node) and 9 p.m. 
(ascending node) local solar time. 


VTPR instruments scan perpendicular to the satellite motion, 
from left to right while facing the direction of travel, in 23 
discrete steps. Fields of view, scanning times, and apparent 
motion on the earth provide "spots" which are contiguous both 
across and along the orbital track. In figure 1 the shaded areas 
indicate the areal coverage during two orbits. Darker shading 
indicates areas where data are redundant during the two orbits. 


Coverage of the earth by VIPR is not complete equatorward of 
49° because the instrument scans only + 30.3° from the local nadir. 
At low latitudes as little as two-thirds of the area is observed, 
but within a 24-hour period the combined ascending and descending 
portions of the orbits provide nearly complete coverage. 


A VTPR instrument observes each spot in eight spectral inter- 
vals 4 cm7! to 16 cm-! (.09 um to .56 um) wide. One interval is 
in the window at 12 um, six are located in the 15 wm band of car- 
bon dioxide, and one is the 19 wm region where water vapor absorp- 
tion dominates. After a set of eight spectral measurements has 
been obtained, the scan mirror is stepped to the next spot. When 
the 23 spots have been observed, the mirror returns to the origi- 
nal position. The pattern resulting from the scanning is depic- 
ted in figure 2, which is an enlargement of the small outlined 
area in figure 1. Observations within this area are obtained 
in about five minutes. 


Analysis of the data first requires their reduction to "clear 
radiances!’ which are spectral radiances that one would find with 
completely clear skies. As described in a later section, data 
are Separated into boxes of adjacent spots, from which a statis- 
tical analysis produces a single group of eight radiances for each 
box. Results are used to retrieve atmospheric profiles. 


The scheme initially applied to the VTPR data is shown in fig- 
ure 2. Eight scan lines are divided into three boxes of 8 x 8, 
8x 7 and 8 x 8 spots. Clear radiances are deduced for three 
boxes indicated by X's in the figure. A retrieval is made for 
each box. The dimensions of the boxes are 500 by 600 km and 700 
by 600 km for the center and side boxes, respectively. New pro- 
cedures will be instituted to locate retrievals at the centroids 
of clear areas. 


Although there are potentially about 200 soundings per orbit 
(2,600 daily), several factors limit the number that can be pro- 
vided to users. These include excessive cloudiness, which makes 
it impossible to obtain satisfactory soundings in the troposphere; 
multi-layered clouds, which reduce the reliability of clear radi- 
ances; high and varied terrain, which produces a similar effect; 
and redundance over polar regions. 


The National Environmental Satellite Service will continue to 
improve the scope and quality of soundings through an active de- 
velopment program. Occasional changes in the methods of deduc- 
ing soundings will be made, and ancillary data, such as those 
from the Scanning Radiometer and from statistical regressions, 
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will be incorporated. As significant alterations are made, and 
aS parameters of newly launched instruments are needed, supple- 
ments to this report will be issued. 


Parts of this report were prepared by the following personnel 
of the National Environmental Satellite Service. 


Section 1 - Introduction (D.Q.Wark) 
Section 2 - Spacecraft (J.M. Siomkajlo and P.G. Abel) 
and Ground Equipment (L.M. McMillin) 
Section 3 - Data Reduction and Analysis 
A. Overall Data Flow (A. Werbowetzki ) 
B. Geographical Location of Data 
(L.A. Lauritson) 
C. Calibration Procedure (J.A. Pritchard, 
D.S. Crosby, D.Q. Wark and A. Werbowetzki ) 
D. First Guess Fields for VTPR Data Proces- 
sing (H.M. Woolf) 
E. Procedure for Obtaining Clear Radiances 
(R. Luebbe, L. A. Lauritson, and 
L.M. McMillin) 
F. Retrieval (M.P. Weinreb and H.E. Fleming) 
Section 4 - Quality Tests and Rejection Procedures 
(F.E. Bittner and C.M. Hayden) 
Section 5 = Data Outputs and Archives (L.M. McMillin) 


2. SPACECRAFT AND GROUND EQUIPMENT 
A. Spacecraft 


The Improved TIROS Operational Satellite (ITOS) is a 1.02 x 
1.02 x 1.25 m, 336 kg spacecraft designed to operate in a sun- 
synchronous polar orbit at 1464 km. Power is supplied by three 
1.65 x .91 m solar panels, which are extended in space. The 
Digital Data Processor (DDP) receives data from the VTPR sensor 
and forwards them to the magnetic tape recorders. A flywheel is 
used to restrict satellite attitude errors about the pitch axis 
to + 0.59, and electrical coils are used to correct roll and yaw 
errors (Schwalb 1972). 


B. Vertical Temperature Profile Radiometer (VTPR) 
(1) Description of the Instrument 
The VITPR is a continuous day-night eight-channel sounding sys- 


tem. Six of the eight channels are used to deduce radiances in 
the 15 um carbon dioxide band. Two other channels are used to 
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deduce radiances at 535 cm-l (in a water vapor absorption band) 
and at 835 cm-l (in an atmospheric window region); these two 
radiances are used respectively to determine the atmospheric 
humidity and to act as a control in the determination of clear 
radiances for the other seven channels (see Section 3.B). The 
system is designed to permit calculation of the vertical temper- 
ature profile from the earth's surface to about 30 km, and to 
obtain an estimate of the moisture in the lower troposphere. 


The optical system of the VTPR instrument consists of a scan- 
ning mirror, a 73.5 mm Cassegrainian telescope, a chopper, and a 
filter wheel (all shown in fig. 3); and a detector assembly 
(not shown) which consists of an Irtran 4 lens, a wide band anti- 
reflection-coated germanium window, and an uncooled pyroelectric 
detector. Broad areal coverage is obtained by scanning + 30.39 
(center of the center scan spot to the center of the outer scan 
spot) from the nadir direction in 23 incremental steps of approx- 
imately 2.7° (fig. 4). At each scan spot, measurements are ob- 
tained for all eight channels as the filter wheel completes one 
revolution. In order to provide image motion compensation, so 
that the measurements for all eight channels cover the same geo- 
graphical area, the filters on the filter wheel are arranged in 
aispival (fig,. 5). 


A black temperature-controlled chopper provides a reference. 
Radiation from the alternating signal resulting from successive 
views of the chopper and the scene is passed through one of the 
filters on the filter wheel and through an Irtran 4 lens to the 
detector. A mask between the lens and the detector limits the 
field of view to a 2.1369 by 2.2369 rectangular scan spot. The 
projection of the scan spot on the earth's surface is approxi- 
mately a square 55 km on a side when the satellite is viewing in 
the nadir direction. Scan spot size increases at larger viewing 
angles because of the increased distance from the earth's surface 
to the satellite and the curvature of the earth (fig. 4). 


Measurements of spectral radiance must be accurate if they are 
to be used to obtain atmospheric temperature profiles. The VTPR 
instrument provides calibrated radiance measurements ranging 
from zero to 204.8 mWXm2 sr cm-l), enabling the instrument to 
measure equivalent source temperatures up to 340 K. Maximum 
allowable relative error between any two channels other than the 
nominal 668.5 cm-l channel is 0.25 mW<m2 sr cm-l). Maximum allow- 
able relative error between the nominal 668.5 cm-l channel and 
any other channel is 0.75 mW m2 sr cm-l) or less. In order 
to avoid certain problems which can occur when analog data are 
digitized on the spacecraft. One digital count corresponds to 
approximately 0.3 mV / (m2 sr cm-l), 
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Analog output from the detector is digitized into 10 bits of 
binary data. Scan spot location is identified by a five-bit 
binary code, which is combined with the radiometric data and a 
parity bit to form a 16-bit word. The resulting digitized data 
from the VTPR are stored in an output register, transferred to the 
spacecraft Digital Data Processor, and stored on tape. Digital 
data from the VTPR and other spacecraft components are trans- 
mitted on one of four channels of the S-band link to the ground. 


A single scan of the VTPR instrument takes 12.5 s, of which 
11.5 s are used to sample the 23 scan spots. Measurements for all 
eight channels are obtained during the 0.5-s interval allowed for 
each scan spot. One second is required for the mirror to return 
to its original position; during this time measurements of various 
components needed for calibration are recorded in place of the 
16 radiance measurements. Output data for a normal scan line are 
summarized in table 1 of appendix I. When the instrument is ope- 
rated in its automatic calibration mode, a 37.5=-s calibration se- 
quence is recorded (replacing 3 scans of normal data) every seven 
minutes. Data from this calibration sequence are used in the cal- 
ibration procedure described in Section 3 C. During this sequence, 
measurements of electrical voltage calibration, space calibration, 
housing calibration, and housekeeping telemetry are digitized and 
inserted into the data stream. During regular operations, the 
calibration sequence is initiated only once per orbit by direct 
command from the ground. 


Table 1 shows the nominal characteristics of the eight filters 
in the VTPR instruments. In column 1, filters 1-6 are numbered 
according to the opaqueness of the atmosphere in the respective 
spectral intervals, so increasing numbers correspond to decreas- 
ing atmospheric opaqueness. This means that the lower numbered 
filters pass radiation containing information about the upper 
atmosphere and the higher numbered filters are used to measure 
the atmospheric temperature near the surface. The numbering sys- 
tem shown in column 1 is used in this report. Numbers in column 
2 show the order in which the filters are mounted on the filter 
wheel. This sequence was chosen to minimize changes in signal 
level between adjacent filters, thereby minimizing the effects of 
any memory errors present in the instrument. 


A summary Of the VTPR parameters is given in table 2. 


Table 1.--Nominal spectral intervals for VTPR filters 


Filter 

Filter wheel 
number position 

1 1 

2 8 

3 2 

4 7 

5 3 

6 6 

7 4 

8 5 


Table 2.--Summary 


Parameter 


Spectral range 

Line rate 

Field of view 
Dynamic range 
Sensitivity” 

Digital signal output 
Data rate 

Primary f/number 
Effective £/number 
Spectral resolution 


Scan-mirror aperture 


Filter band pass 


Center wavelength Half-width Tenth-width 
(wm) 


14. 


14. 


14. 


14. 


fig 


LoS 


18. 


ea 


96 


aM) 


Go 


o7 


(cm71) (cm-l) (cm-1) 
668.5 35 TOSS 
677.2 10 20 
695.0 10 20 
708.0 10 20 
I2o20 10 20 
747.0 10 20 
535.0 18 20 
$33.0 10 16 


of VTPR parameters 


Nominal value 


12-19 um 

4.8 lines /min 

2.1362: by. 2.236" 

0 to 210 mW/(m* sr cm!) 
0.25 mW / (m2 sr cem7+) or less 
16 bits 

256 bits /s 

£73 

£/0.6 at the detector 

See table 1 


3 in. 


*O~Branch 0.75 mV / (m2 sx cm-l) 
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(2) Spectral Transmission Functions of VTPR Filters 


All VTPR filters, with the exception of the Q-branch filters at 
668.5 cm-l, were manufactured by the Optical Coating Laboratory, 
Ine, (OCLI). Q-branch filters were made by Grubb-Parsons, Ltd. 
(GPL), and’ by Spectrum Systems, Inc. (SSI). Table 2 of appendix 
I gives the serial numbers assigned to the VIPR instruments on 
NOAA 2, as well as the corresponding best estimates of filter 
spectral characterisitcs. 


The importance of obtaining accurate spectral data and the 
difficulty of doing so implies that the manufacturer's original 
data should be verified whenever possible. Data for all instru- 
ment 1 filters, except the nominal 677.5 cem71 filter, were ob- 
tained at the NESS Satellite Experiment Laboratory (SEL). This 
was, unfortunately, not possible for those in other instruments, 
for which data on only the 535 cm-l filters were obtained at SEL. 


In the SEL verification procedure each of the VIPR filters was 
mounted in a special holder at 7.5° away from the normal to the 
£/10 beam in a Beckman IR-7 spectrometer. The_holder was rigid- 
ly mounted to the spectrometer so that the 7.5° angle was accu- 
rate to better than + 0.39. Filter temperature was set near 
35C, typically within + 0.5C, and the instrument was purged con- 
tinuously with pure dry nitrogen. Single beam mode, with a 
spectral resolution of between 0.8 and 1.0 cm7~!, was chosen for 
all the runs. This resolution causes a small but significant 
convolution distortion in the case of the Q-branch filter only. 
However, this did not seriously affect the comparison with the 
manufacturer's data, which were similarly distorted. The distor- 
tion is always in the direction of a less desirable filter pro- 
file. Thus, filters whose measured profiles meet the specifica- 
tions are acceptable. 


Spectral calibration of the spectrometer was achieved using 
known carbon dioxide absorption line positions between 577 and 
959 cm-l. For the latter region of weak absorption, a White cell 
with an optical path length of 10 m was used. The calibration 
runs were made as closely as possible in time to the data runs, 
but sometimes there was an overnight delay. 


The spectrometer output was recorded in analog form on a chart 
recorder and in digital form on a tape recorder. The tape record 
was processed to provide a computer-—~drawn plot of transmission as 
a function of indicated (uncalibrated) wavenumber. 
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Table 3 of appendix I presents comparative figures on spectral 
Characteristics obtained from all sources of data. Agreement is 
always within the noise level of the measurements except for the 
535 cm-l filter. In this case, there is a significant disagree- 
ment in center wavenumber, most probably caused by an unadvoid- 
ably incomplete calibration of the SEL spectrometer. The OCLI 
figures have therefore been accepted for the center wavenumber of 
that filter. In addition to using CO» line positions near 690 
cm7!, the OCLI spectrometer was also calibrated with the water 
rotation line at 525.97 cm-l! (vacuum). 


Appendix II contains the best estimates of the filter spec- 
tral characteristics in instruments 1 through 4. They are uncor- 
rected for convolution distortion. The curves were calculated 
solely from the OCLI chart traces. Points were taken from the 
curves about every 0.4 cm! and smooth curves were fitted to the 
points by means of spline functions. 


C. Other Spacecraft Equipment 


The NOAA 2 satellite carries a two-channel Scanning Radiometer 
(SR), a Very High Resolution two-channel scanning Radiometer 
(VHRR), and a Solar Proton Monitor (SPM) in addition to the VTPR 
instrument (Schwalb 1972). Sea-surface temperatures, which are 
used in the VTPR data processing, are derived from SR measure- 
ments covering the 10.4-12.5 um spectral interval. The SR in- 
strument views a smaller area than the VTPR instrument and thus 


has a better chance of receiving radiation from a cloud-free area. . 


Statistical techniques are used to identify the measurements that 
are cloud-free, and the cloud-free radiances are used to deduce 
the ocean surface temperature. 


D. Ground Equipment 


Ground equipment required for the VTPR data processing consists 
of two readout stations (Wallops, Virginia and Gilmore, Alaska); 
communication lines connecting the readout stations to the NOAA/ 
NESS data processing and analysis facility in Suitland, Maryland; 
and the following types of processing equipment and computers at 
Suielands 


RCA Modulator - Demodulator (MODEM) 2270678-501 
EMR 2721 -— signal conditioner 

EMR 2731 - PCM frame synchronizer 

EMR 6130 — computer 

CDC 6600 = computer. 


LS 
3. Data Reduction and Analysis 
A. Overall Data Flow 


The VTPR data flow and software system consists of six major 
computer programs. Figure 6 shows the individual program modules 
and their relationship within the flow system from the spacecraft 
to data users. 


While the NOAA 2 satellite is above the horizon of either of 
the two readout stations, the digital data are played back from 
an onboard Scanning Radiometer Recorder (SRR) and immediately for- 
warded via communication lines to the NOAA/NESS data processing 
and analysis facility in Suitland, where they are demodulated by 
the MODEM and passed on to the EMR computers for initial proces- 
Sing. The VTPR data are separated from the other spacecraft data 
and finally written on an output digital magnetic tape by the 
Ingest Program in the EMR 6130 computer. Normally, one orbit of 
data is processed for each readout. However, once each day the 
satellite does not appear above the horizon at either of the two 
readout stations during one or two orbits, and it is necessary to 
store these orbits of data on the spacecraft until the next read- 
out. In these cases, up to 209 minutes of VTPR data may be re- 
corded on one output tape, and two tapes may be required to re- 
cord all the data. 


A sequence of programs is executed on each output tape from the 
EMR 6130. The first program computes the geographical location 
of each vertical sounding to be produced. Data concerning the 
instrument status and condition are then analyzed, and parameters 
such as electrical gain and standard deviations are checked and 
verified to be within acceptable limits. Following the proce- 
dures described in section 3.B, radiance values are computed and 
passed, together with appropriate earth location and time infor- 
mation, to the next program. An archive tape containing the 
earth-located radiances also is generated. 


The Clear Radiance Program operates next to eliminate the 
effects of clouds. It statistically computes an equivalent clear 
radiance for an area by comparisons of adjacent scan spots, uti- 
lizing knowledge of sea-surface temperature (obtained independent- 
ly from the Scanning Radiometer data) and an approximate first 
guess vertical temperature structure (obtained from climatology 
and a recent forecast or analysis). Earth-located clear radi- 
ances are generated and added to the archive tape. The data are 
then used to retrieve temperature profiles. 
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The Retrieval Program uses the first guess temperature profile, 
a first guess humidity value, and atmospheric transmittance val- 
ues to compute the vertical atmospheric profiles from the input 
of clear radiances. 


The final output data are quality checked and reformatted onto 
three magnetic tapes. One tape is forwarded to the National 
Meteorological Center (NMC) for input to their numerical analysis 
and primitive equation (PE) forecast model. The second tape is 
given to the National Weather Service communications center for 
transmission to users. Data that pass the quality check are 
added to the archive tape, which is then transmitted to other 
users (Goddard Institute for Space Studies, for example) and to 
the permanent archive of the NOAA Environmental Data Service. 


B. Geographical Location of Data 


Earth location of the data is performed in two steps. The first 
step is performed in an Earth Location Program, which calculates 
latitudes and longitudes for centers of the areas for which VTPR 
temperature and humidity profiles are to be calculated. fThe sec- 
ond step is performed as a subroutine of the calibration program. 
This subroutine calculates the locations of the VIPR scan spots 
in the grid on which the sea-surface temperatures are provided. 


The location of a scan spot is determined by a point called the 
sensor principal point. This point is the center of a scan spot 
on a plane surface perpendicular to the viewing angle of the in- 
strument. Because of the curvature of the earth, this point is 
the exact center of a scan spot only when the instrument is view- 
ing along the local nadir. In the earth location procedure, prin- 
Cipal points are calculated for spots 5, 12, and 19 of the fourth 
of eight scan lines (fig. 9). To locate these principal points, 
the satellite position, the satellite attitude, the sensor mount- 
ing position, and the mirror position must be known. Attitude 
information in the form of max and A (RCA 1972) is supplied 
by the Satellite Operations Control Center (SOCC); $a, is the 
maximum roll angle during an orbit; and ’} is the angle from the 
ascending node to the point of maximum yaw. Values of $5, and 
A} are easily converted to roll and yaw angles. The SOCC attempts 
to control &)4x to within 0.5° tolerance, with a maximum error of 
about 1°. The pitch tolerance is generally kept to within 0.5°. 


Weekly orbital element data are supplied by the NASA. These 
Gata and the General Electric orbital predictor package (Brower 
1959, and Lyddane (1963) are used to determine satellite position. 
Once the satellite position is determined, a series of matrix 
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rotation and vector additions are used to compute the principal 
point latitude and longitude (figs. 7 and 8). These rotations 
convert the geocentric principal point to a position in a geodet- 
ic coordinate system, and include corrections for satellite alti- 
tude, satellite position, sensor mounting position, and mirror 
position. 


Sea-surface temperatures from the SR are available for i and j 
points of a square grid superimposed upon a polar stereographic 
projection. Positions of the bench marks (fig. 9) are calculated 
for the VTPR location and converted to positions on the polar 
grid. Temperatures of intermediate boxes are obtained from the 
polar grid using coordinates obtained by interpolation from the 
known bench mark coordinates. 


C. Calibration Procedure 


Calibration of the VTPR sensor data involves two procedures. 
The first is a pre-launch ground calibration of the instrument; 
the second is used after launch to adjust the calibration for 
Changes in the spacecraft operating conditions. 


During the first procedure each VTPR instrument is calibrated 
by being exposed to an external reference source. The internal 
source iS not a standard and must therefore be calibrated against 
a standard source. This is accomplished by calibrating the in- 
ternal source with an external source which, in turn, has been 
Checked with. a standard. 


Calibrations of the VTPR are based on the external sources used 
during thermal-vacuum tests on the ground. The honeycombed ex- 
ternal sources used in the calibration have provisions for main- 
taining any fixed temperature between 180 K and 340 K. Six ther- 


mocouples are used to measure the temperature; they indicate that 


thermal gradients are less than 0.1 K across the surface of the 
source. These external sources have been compared with the low 
temperature blackbody reference of the Canadian National Research 
Council (NRC) (Bedford 1970), using calibrated thermopile detec- 
tors as transfer standards (Hilleary et al. 1969). Thermopile 
output for the VIPR references read 1 + 0.5% higher than the out- 
put from the NRC reference. 


During the calibration of the VTPR internal source, two of the 
external sources were placed in the "earth" position so that they 
filled the field of view for scan spots 1 and 23, the extreme 
mirror positions. With the instrument temperature and the scan 
position fixed, a calibration was performed by varying the 
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A = ascending node of orbital plane 

C = north pole 

F = geocentric subsatellite point 

I = satellite position 

AFE = geocentric subsatellite point track 

JHFC = meridian through F 

BOD = right ascension of orbital plane 

COE = inclination of orbital plane 

AOF = Q(angle from ascending node to geocentric satellite 
point 

FOH = 82 and AFH = 86), angles to correct geocentric subsatellite 

point to geodetic subsatellite point 
FOJ = geocentric subsatellite latitude 
HGJ = geodetic subsatellite latitude 


Figure 7.--Location of satellite position 


HIK = roll, pitch, and yaw rotation 

K = principal point of a VTPR data spot 

KOM = geocentric latitude of principal point 
KLM = geodetic latitude of principal point 
CKM = mderdian through K 


Figure 8.-- Location of principal point showing effect of 
roll, pitch, and yaw error 
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temperature of one of the external sources in steps between 180 K 
and 340 K. This procedure was repeated for a series of instrument 
temperatures. During this calibration process, the VTIPR was put 
into its automatic calibration mode a number of times to allow 
comparison of the internal source with the external secondary 
standards. 


The second procedure is performed when an instrument is in 
space. Then, the only radiation references are the internal 
blackbody source at the nominal temperature of the instrument 
and a view of space. At the wavelengths used in the VTPR, 

a space view effectively has zero radiance. For a linear system, 
two references provide sufficient information to calibrate the 
output of the instrument in physical units when it is viewing the 
earth. A relationship between source radiance and the digital 
output of VTPR is given by 


4 
Ti; = aj + by Ci (1) 


= radiance 

= calibration coefficients 

VTPR output in counts 

spectral interval 

= source identifier (space or internal), 


Ub Q OH 
Il 


The linearity of the output implied by eq (1) is an integral 
part of the instrument design, and has been verified by tests of 
each instrument. 


In the normal earth-directed views, spots 1-23, the coefficients: 
"a" and "b" for the ith channel are determined from 


3 
Se 2 ae ae T 
if K=0 ik k t 
and A 
2 1 
ge = ee a 
1 K=0 ik k 
where T. = 1.0 t 
0 
T) = primary optics (counts) ; 
T, = secondary optics (counts) t 
q 
Tz = shroud (counts). 


The last three terms in eq (2) account for the contributions of 
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components of the VTPR to the radiative flux at the detector. 
Emissions of these components are temperature dependent, but be- 
cause of the narrow range of operating temperature for these com- 
ponents, their contributions to the flux may be taken to be a 
linear function of temperature. Values of the coefficients Aj, 
and Biz are determined by regression from data obtained during 
pre-flight thermal vacuum tests, using output counts from the VTPR 
and the various thermal probes listed in eq (2). In space, values 
of T, for these components are obtained for each scan during the 
retrace period. Since these temperatures are obtained during each 
scan, the values of aj and bj can change during an orbit to com- 
pensate for changing conditions on the satellite. 


The values of aj and bj; given by eq (2) are periodically checked 
while the instrument is in space. During calibration checks in 
space, the VTPR views space, identified as spot 27, and the inter- 
nal source ("patch"), spot 28. These views provide data for the 
relations: 


1 a de 27 
and (3) 
Bane 0 85 ee Oe oe 
my = di . 
where B(v, T 31) Planck radiance at v5 and Toy 

vj; = spectral interval 
T 31 = temperature of the "patch" 
fal Tal = : "W "W 
Cc. 97 and Cc; 23 7 average counts during the "space" and 


"patch" views at v;. 


This procedure provides a check of the values of aj; and b,. When 
the values do not check the coefficients in eq (2) are adjusted. 
As a last resort the instrument may be placed in an automatic cal- 
ibration mode, values of aj and b; will be determined from eq (3), 
and new coefficients will be generated for eq (2). 


Data for the calibration check are provided during a calibra- 
tion sequence. This sequence (table 4 of appendix I) lasts for 
37.5 seconds. In the automatic calibration mode, the instrument 
obtains calibration data approximately every 7 minutes throughout 
the orbit. Normally, the VTPR is commanded to obtain calibration 
data only once per orbit. During a calibration sequence, an 
electrical calibration is performed using the staircase voltage 
counts to verify the linearity and stability of the electronic 
circuitry. Also, means and standard deviations are computed for 
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each phase of the calibration sequence (electrical staircase, 
space view, and patch view). 


Instrument temperature data recorded during the one-Second re- 
trace period are identified by a five-bit binary code. Prior to 
launch, each instrument was arbitrarily assigned identifier code 
24 or 25. During the one-second retrace interval of each mirror 
scan, Output counts from the probes measuring the primary optics, 


secondary optics, and shroud temperatures, a 24-bit time code, and 


a frame synchronization identifier are recorded. 
D. First Guess Fields for VTPR Data Processing 


When describing the nature and source of first guess fields 
utilized in VTPR data processing, it iS convenient to divide the 
earth into three geographical regions: (1) that portion of the 
Northern Hemisphere north of 18°N (hereafter identified as NH); 
(2) the area between 18°N and 18°S (hereafter identified as 
tropics); and (3) the Southern Hemisphere south of 18°S (SH). 


(1) Northern Hemisphere 


Guess fields for the NH are extracted from the operational fore- 
cast of the National Meteorological Center. The fields used and 
the time intervals at which forecasts are available are listed 
in table 5 of appendix I. 


When forecasts are available for a given field, the forecast 
nearest in time to that of the satellite observation is used. 


The fields are read from system disk (permanent files) and 
stored in Extended Core Storage (ECS). To construct a guess pro- 


file for a particular VTPR observation point, the fields are read | 


into central memory from ECS, and a value for each parameter is 
interpolated to the VTPR location from the four surrounding grid 
points. The portion of the guess temperature profile above 10 mb 
is generated by regression as explained in 3. F. 


(2) Tropics 
For the tropics, radiosonde data averaged over an appropriate 
two-month period provides the guess temperature profile; regres- 


sion estimates are employed to estimate the relative humidity. 


(3) Southern Hemisphere 


In the Southern Hemisphere, a feedback technique is used to gen- 


erate guess temperature fields up to 10 mb. The retrieved VTIPR 


25 


profiles for a given day, together with other available upper air 
data, are used as input to a multi-level objective analysis pro- 
gram, whose output,in the form of grid point data, provides the 
first guess fields for the subsequent day's VTPR data processing. 
Humidity information is provided by regression estimates. 


E. Procedure for Obtaining Clear Radiances (Program CLRAD) 
(1) Mathematical Justification 


To eliminate the effect of clouds on the retrieved temperature 
profiles, a single cloud-free or "clear" radiance is produced 
from measurements obtained from a number of scan spots. The re- 
sulting "clear" radiances are then used to obtain temperature pro- 
files. "Clear" radiances for a 7 x 8or an 8 x 8 box or sSubarray 
are obtained from the measured radiances of the scan spots and the 
sea-surface temperatures provided by data from the Scanning Radi- 
ometer (SR). They are obtained using a technique suggested by 
Smith in Fritz et al. (1972), and in Smith et al. (1970). 


For an area that is partly cloud covered, the radiance can be 
expressed as 


I(v,) = NIK(v;) + (1-N) T1 (¥y) (4) 


where N is the fractional cloud cover, Icy (vi) is the radiance of 
the cloudy areas, I,j],(vj) is the radiance of the clear areas, 

and v; is the central wavenumber of spectral interval i. The term 
Toy (vi) Can be written as 


Dey i? = € (vi Tg (%y) + [1l-e(v,)] Tar (Yy? (5) 


where €(v;) is the cloud emissivity and Tceld(vj) is the radiance 
of a completely opaque cloud. Using eq (5) to substitute for 
Toy (vz) in eq (4) leads to 


I(v,) = a(v. Ig (yy) + [l-a(v. ) IT (yy) (6) 


where @(v;), the product of the fractional cloud cover and the 


emissivity, is given by 


= Z 
a(v,) = Ne (v,). (7) 


In the spectral region covered by the VTPR, clouds can be assumed 
to be grey, so a is therefore independent of wavenumber. If there 
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are two adjacent scan areas with different values of @ but with 
the same values of Igjq(vi) and Ieir (vy), then, for a given spec- 
tral interval, a separate equation can be written for each of the 
two adjacent areas 


le @, (vs Tag (3) or [l= (v, JIT, (yy) (8) 


and 


v.) = a, (vi Tg (Ys) + [1-@, (v, IT (yy), (9) 


where the subscripts 1 and 2 denote two areas. Elimination of 
Iciq(vzy) from eqs(8) and (9) gives 


-l 


The atmospheric transmittance for spectral interval 8, which is 
in an atmospheric window, is close to unity. The radiance for 
this spectral interval is only slightly affected by the atmosphere 
SO a very good approximation of the value of Igj;7(vg) can be com- 
puted using the first guess atmospheric profiles of temperature 
and humidity and the surface temperature obtained from the SR data 
With a value of Igir(vg), the ratio a9 (vq) /ey (v3) can be obtained 
from eq (10). Since the ratio is independent of wavenumber, the 
value of the ratio obtained for vg can be used in eg (10) to ob- 
tain values of Icir(v,) for the other spectral intervals. Writ- 
ing eq (10) for vg and some other vj, and eliminating the ratio 
a (vs) fay (vi), gives 


oi¥)=to (99) Tag 9,2 E Oy? 


Tear Mg) Ta (%g) Tong (Vg) - Ty (¥Q) (T1)§ 


The right hand side of eq (11) is the slope of a line containing 
the points [Icir(vjy),Icir(vg)] and [1] (v_),Ig(vg)]- Graphically 
(see fig. 10), the three points [Igjy(¥j),Teoir (vg) 1, [11 (vz), T1 (vg) 
and PCy My T9 (vg) ] lie on a straight line that can be determined 
from the Passuecd values of Ij (vy),11 (vg), I2 (vy) and Ig (vg). The 
value of Igjiy (vj) can then be found from the known value of 

Toelr (va). When calculating values of Icir (vj), radiances from 
adjacent spots with different nadir angles are adjusted to a com- 
mon nadir angle as described later. 
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Figure 10.--Procedure for determining 
clear radiances 


(2) Operational Procedures 


To reduce the noise level of the deduced clear radiances, 
Icir (vi), values from a number of adjacent pairs are combined. 
Figure 2 shows the scanning pattern of the instrument. Data from 
eight successive scans are used to obtain three temperature pro- 
files. Figure 2 also shows how the 23 x 8 scan spots are divided 
to form two 8 x 8 subarrays and one 7 x 8 subarray for which 
three temperature profiles are calculated. 


Scan spots are considered to be adjacent if they meet at a cor- 
ner or an edge. For a given scan spot, pairs are obtained for 
four adjacent scan spots: the one above and to the right, the 
one directly above, the one above and to the left, and the one | 
immediately to the left of the given scan spot (fig. 11). Ineach 
subarray, there are 49 scan spots for which the four combinations | 
with adjacent spots shown in figure 11 are possible. As a re- 
sult, 49 x 4 or 196 values of Igiy(vj) are obtained for each sub- 
array. Except for scan spots on the edges of the 23 x 8 array, 
each scan spot has been compared with each of its eight adjacent 
scan spots when the array is completed. 


The 196 values of Igiy (vj) are combined to obtain an estimate 
of the value of Icir(v;i) for the area. Several procedures for 
obtaining a combined value have been tried using simulated data. 
The best results were obtained by using a mixture of two of these 
procedures. 


When a Single cloud layer was simulated, the resulting distri- 
bution of values of Icir(vj) was symmetrical, and a weighted 
average gave the best estimate of the true value of Tap y (vio. 
The weights used are proportional to the inverse of the estimate 
of the variance of the I,jy(v;) for a given scan spot pair and 
are given by 


M 
Wee, 2 ca (12) 
m m 

ins 


= 
II 


G 
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2 2 
nee [Ty (vg) 25 (vg) 12-7 A Eo, (Vg) Ty (vg) 


2 
#1Z 4, (Yg)-Ty (vg) 14 - (13) 


In eqs (12) and (13), m refers to one of the 196 possible pairs, 
and M is the number of pairs included in the average. 
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Figure 11.—Scanning pattern showing 
the pattern used for adjacent pairs 
and the nominal location of 
retrievals 
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When several cloud layers were used in the simulation, the dis- 
tribution became asymmetrical and the mode gave the best estimate 
of the true value of I,j+(vji).- However, when the mode is deter- 
mined from the number of occurrences in a given class interval, 
an uncertainty occurs because of the size of the class interval. 
To avoid this uncertainty, the distribution of values of Ig1;7(vj) 
was convoluted with a smooth function and the mode was taken at 
the maximum value of the resulting function. Several smoothing ~ 
functions were tried, and a fourth degree chi-square distribution 
was selected. Figure 12 shows typical unsmoothed histograms and 
the smoothed histograms that result from convoluting the data 
with a chi-square distribution. Since the effect of clouds (and 
thus the corrections) are small for channels 1 and 2, the weighted 
averages are always used and histograms are not produced. 


Selection of the value from the weighted average or the mode is 
determined by comparing the two values. Agreement between the 
two indicates that the distribution is symmetrical; in such cases 
the weighted average is selected. Otherwise the mode is selected. 
At the present time, a difference of l mwW/ (m2 sr cm-l) or less 
between the two values is considered to indicate agreement and 
thus a symmetrical distribution. 


During processing of the VIPR data, two checks are made while 
calculating values of Igiy(vi). First, if the measured window 
(vg) radiance "agrees" with the window radiance calculated from 
the first guess, the spot is assumed to be clear. Window radi- 
ances are considered to be in agreement when the measured window 
radiance equals or exceeds the calculated value. Values of Iqir 
(vi) from clear spots are averaged. If one or more clear scan 
spots is found, the value of Igir(vj) is obtained exclusively 
from the "clear" areas. As a second check, values of window ra- 
diance for adjacent scan spots are compared. When values are too 
close, the resulting estimate of Igiy(vj) is considered to be un- 
reliable unless one or more of the spots is clear. Values of 
Iclr (vj) resulting from adjacent scan sage with differences in 
value Of I(vg) that are less than 1.0 mW/(m* sr cm-l) are not 
used. If no clear areas are found and if fewer than 25 values 
of Igiy (vj) are obtained for a 7 x 7 subarray, a retrieval is not 
attempted. 


Measurements for some adjacent pairs are obtained at different 
viewing angles of the earth. When these measurements are com- 
pared, a correction must be made to compensate for the different 
atmospheric path lengths. When comparing radiances in adjacent 
scan spots, and when averaging the results over a7 x 7 subarray 
for adjacent scan spots, the first guess temperature profile and 
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Figure 12.-- Typical histograms and 
smoothed histrograms of radiances 
for channels 3 through 7 
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the sea-surface temperature are used to calculate spectral radi- 
ances for the four cases of no clouds and of 50% cloud cover at 
700, 470, and 300 mb. Radiances are calculated for several zenith 
angles, and curves of radiances verses zenith angle are calculated) 
for the four cases (fig. 13). Calculated and measured first guess 
window radiances are compared to determine which of the four 
curves gives the best estimate of cloud height (fig. 13a). The 
corresponding curve is then used to determine zenith angle correc- 
tions for the other spectral intervals (fig. 13b). When correc- | 
ing to the center of the array, the curve for no clouds is used 
Since clear radiances derived in program CLRAD are being adjusted. 
While the first guess atmosphere is used in making corrections | 
for zenith angles, the zenith angle corrections are small, and 
the first guess atmosphere has a very minor influence upon the 
radiances passed on to the retrieval program. One should also 
note that the value of Loir (vg) resulting from the clear radiance 
procedure is the value calculated from the first guess. 


F. Retrievals 
(1) Introduction 


The temperature and moisture retrieval program derives atmos- 
pheric temperature and humidity profiles from the clear radiances 
produced by the CLRAD program and can be described as a sequence 
of five sections called Input, Construction of Transmittances and 
Weighting Functions, Temperature Retrieval, Moisture Retrieval, 
and Output. 


(2) Input 


The retrieval program uses the radiance values for the first 
seven of the eight spectral intervals listed in table 1. As ex- 
plained in section 3E., the eighth radiance value is calculated 
from the first-guess surface temperature and first-guess atmos- 
pheric profiles. The retrieval program uses the surface tempera- 
ture in place of this radiance. 


A first-guess temperature and water vapor profile is required 
by the retrieval program. Below 10 mb, the first-guess profiles 
are obtained from forecasts, analyses, and climatology as ex- 
plained in section 3 D. The tropopause temperature, height, and 
pressure are also provided by the National Meteorological Center 
(NMC). Above 10 mb, the first-guess temperature profile is ob- 
tained by regression on radiances measured in spectral intervals 
1 and 2 and the first-guess temperatures at 10, 30, and 50 mb. 
This technique gives reliable estimates of the temperature above 
10 mb because it is based upon a recent sample of rocketsonde and 
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Figure 13.--Calculated window radiance versus zenith angle for 
50% cloud cover at four levels showing the zenith angle 


correction 
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radiosonde data (Gelman et al., 1972). It is especially 
valuable in accounting for sudden warming in the stratosphere. 


Finally, the guess temperature profile is interpolated to the 
100 pressure levels shown in table 6 of appendix I. These levels 
are distributed between 1000 mb and 0.01 mb in even increments on 
a scale of pressure raised to the 2/7 power. 


A guess for the water vapor mixing ratio profile is also re- 
quired. In the Northern Hemisphere, NMC provides an average rel- 
ative humidity for the atmosphere between the surface and a vari- 
able level which is near 500 mb. From this single quantity a 
mixing ratio profile is specified as follows: a saturation mix- 


ing ratio profile is calculated from the guess temperature profile. 


between the surface and 500 mb. In this region, the guess mixing 


ratio profile is constructed as the product of the relative humid-_ 


ity and the saturation mixing ratio. Above the 500 mb level, the 
mixing ratio profile is constructed to fall off as p>. 


In areas where no moisture information is provided by NMC, such 
as the Southern Hemisphere, the correlation between temperatures 
and mixing ratio is utilized through the regression relation 


we -W= 4H (T* - T), (1 4) 


where W* is the required guess mixing ratio profile and T* is 
the guess temperature profile. The quantities T, W, and H are 
the mean temperature profile, the mean mixing ratio profile, and 
a matrix of regression coefficients, respectively, which have 
been compiled in advance from conventional soundings. 


(3) Construction of Transmittances and Weighting Functions 


Atmospheric transmittances for the eight spectral intervals of 
table 1 are fundamental to the retrieval and must be constructed 
accurately. The total transmittance of the atmosphere is assumed 
to be the product of the individual transmittances of carbon di- 
oxide, ozone, and water vapor, 


T(p) = (po) een to) re T 9 (P): (15) 


Tco 
2 3 2 


Transmittances for COj are based upon calculations made with the 
point-by-point method of Drayson (1971). His method was used to 
calculate transmittances for the U. S. Standard Atmosphere, 1962 
(COESA 1962) and for atmospheres which differed from the Standard 
Atmosphere by + 10K, + 20K, and + 30K. Temperature corrections 
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to the transmittances, generated from these calculations, are 
used in the retrieval program to adjust the transmittances to val- 
ues appropriate for the first-guess temperature profile. 


Appendix II contains the measured filter characteristics of the 
first four instruments. It includes carbon dioxide transmittances 
for a typical atmosphere and weighting functions (derivatives of 
the transmittances with respect to p2/7), which indicate the rel- 
ative contributions of levels of the atmosphere to the radiance 
received at the satellite. Radiation detected by the VTPR origi- 
nates from layers which increase in pressure as the channel wave- 
number increases. 


Ozone transmittances are a minor correction and are applied as 
a Single correction for all latitudes and seasons. The transmit- 
tances were calculated for an average profile (McClatchey 1972) 
from a line-by-line technique for pressures greater than 100 mb, 
and from a band model for lower pressures. 


Water vapor has a marked influence upon atmospheric transmit- 
tances in the lower troposphere. The infrared spectrum of water 
vapor includes both spectral lines and a strong continuum, which 
are treated Separately. The line contribution cannot be calcula- 
ted by the same procedure used for CO» transmittances because the 
mixing ratio of water vapor is highly variable. Instead, the 
first step-(performed in advance of the retrieval) is to fit a 
model to line by line calculations of spectral line absorption 
averaged over the appropriate filter functions. These calcula- 
tions are performed over homogeneous paths and cover, in incre- 
ments, the expected atmospheric range of pressure, temperature, 
and mixing ratio. In the retrieval, the atmosphere is treated as 
a succession of homogeneous layers (Weinreb and Neuendorffer 1973). 


The absorption by the water vapor continuum is observed to vary 
as the square of the water-vapor concentration and as the inverse 
fifth or sixth power of temperature (Burch 1970 and Bignell 1970). 
This behavior, which is consistent with absorption by the water 
vapor dimer (H20)9, is included in the continuum transmittance. 

A term ascribed to foreign-gas broadening of distant lines, which 
is proportional to the first power of the moisture concentration, 
is also included. The overall water-vapor transmittance is 
finally calculated as the product of the line and continuum 
transmittances. 


(4) Temperature Retrieval 


The temperature profile is obtained through a modification of 
the minimum-rms solution (Rodgers 1970 and Strand and Westwater 
1968) of the radiative transfer equation. This solution can be 
written as 


BL. (T) = BL (T*) + (C) (R-R_*) (16) 4% 


where the matrix C is given by 


c= sat (asa + n)7, (17)9 


A description of all terms in these equations follows. Eq (16) 
and (17) yield B;(T), which is the 100-element Planck-radiance 
profile of the solution temperature profile, T, computed at a 
reference wavenumber of 700 cm-!. The 100 elements correspond to 
the 100 atmospheric pressure levels shown in table 6 of appendix 


I. A solution temperature profile is readily obtained from B,(T)./ 


Geopotential thicknesses relative to 1000 mb are computed from 
the perfect-gas and hydrostatic equations. In the Northern Hem- 
isphere, heights are obtained relative to the NMC forecast 850-mb 
height. 


The term B,;(T*) is the Planck-function profile computed at 700 
cm-l from the guess temperature profile T*. The vector R, con- 
tains six elements, the measured radiances in the first six spec- 
tral intervals of table 1. Likewise, R*# is a vector of six radi- 
ances calculated (by the radiative transfer equation) from the 


surface temperature, the guess temperature profile, and the trans-. 


mittances. Both Ry and R¥ are scaled to the 700 cm-l reference 
wavenumber through radiance-equivalent temperatures. The 6 x 100- 
dimensional matrix A consists of the six 100-dimensional weight~ 
ing functions illustrated in appendix II. 


Uncertainties in measurements of the radiances, originating 
from instrumental noise and calibration errors, quadrature errors, 
and uncertainties introduced in obtaining clear radiances, reside 
in N, a 6 x 6 dimensional variance-covariance matrix. Statistics 
of the atmosphere comprise the 100-dimensional diagonal matrix S, 
whose elements are variances of the 700 cm-! Planck-radiance pro- 
file derived from a set of typical radiosonde measurements. The 
off-diagonal elements of S, which should contain the correlations 
among the temperatures at different levels of the atmosphere, are 
ignored for computational speed and simplicity. 


a bel Pet en eee eee 
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Eq (16) and (17) express the solution as the guess profile plus 
a linear combination of differences between observed radiances 
and the radiances calculated from the guess profile. Large errors 
in the measurements of radiances (large N) will force the solu- 
tion toward the guess. The solution will also tend toward the 
guess for layers of the atmosphere with small expected variations 
in the temperature (small S). Usually one application of eq (16) 
and (17) will provide a convergent solution in the sense that 
radiances calculated from the retrieved profile differ from the 
measured radiances by less than the standard deviation of the 
errors. Otherwise, eq (16) and (17) are applied again, with T* 
and R# representing the result of the previous iteration. Since 
the matrix C is nearly independent of temperature, it is held 
constant from iteration to iteration, thus minimizing computation 
time. 


Figures 14 and 15 compare operational VTPR retrievals with near- 
ly simultaneous radiosondes for a high latitude (fig. 14) and a 
tropical (fig. 15) case. It is seen that the VTPR retrievals 
follow the trend of the radiosondes up to and beyond 50 mb. How- 
ever, because the weighting functions are broad, the solution does 
not contain fine-scale structure unless it is introduced via the 
guess profile. 


(5) Moisture Retrieval 


The moisture retrieval begins with a measured radiance at 535 
cm-! and the retrieved temperature profile. The solution for the 
mixing ratio profile W is assumed to have the form 


W = W* + Co (18) 


where W* is the guess mixing ratio profile, C is a constant to be 
determined, and » is an empirical orthogonal function (Alishouse 
et al. 1967) computed in advance of the retrieval from conven- 
tional sounding data. If it is impossible to compute 9» because 
of lack of reliable data over certain areas, » is assumed to be 
W* itself. The constant C is evaluated in the retrieval through 
the requirement that 535 cm-! radiance computed from W be equal 
to the measured radiance within two standard deviations of the 
measurement errors. Finally the solution W is transformed to dew- 
point depression. Only one spectral interval is used to measure 
water vapor. With this single measurement, only one parameter 
(e.g., the value of C) can be determined. The distribution of 
the water vapor is determined by the relatively smooth functions 
W* and ». Another limitation of any moisture retrieval is the 
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Figure 14.-- Comparison between a VTPR sounding and a radiosonde 
at 52.59 north, 20° west. D.P. indicates dewpoint. 
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| Figure 15.--Comparison between a VTPR sounding and a radiosonde 
at San Juan (18° north). D.P. indicates dewpoint. 
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dependence on accurate knowledge of the temperature profile 
Since temperature enters strongly into the computation of radi- 
ance. Operational moisture retrievals appear in figures 14 and 
15; the lack of detailed structure in the moisture retrievals is 
apparent. 


(6) Output 


For each sounding, the following quantities are passed to the 
quality control program: 


a. Retrieved temperatures, heights, and guess temperatures 
at the 15 standard levels (1000, 850, 700, 500, 400, 
300, 250, 200, 150, 100, 70, 50, 30, 20, and 10 mb) and 
at the tropopause 


b. Clear radiances in the eight spectral intervals 
(table 1) 


c. Dewpoint depression at the five standard levels up to 
the 400 mb level 


d. Two coefficients to account for atmospheric attenuation 
(limb darkening) in the Scanning Radiometer measurements 
of equivalent temperature which are used to derive 
surface temperatures. dod 


4. QUALITY TESTS AND REJECTION PROCEDURES 


The final output from the VTPR is called the Satellite Infrared 
Sounding (SIRS); it includes location, heights, temperatures and 
dewpoint-temperature depressions at the mandatory pressure levels 
up to 400 mb, height and temperatures up to 10 mb, and tempera- 
tures at significant levels. Initially, only data from open 
water oceanic regions are being released to the user. Before the 
data are released, several tests are performed. Failure of any 
one test causes all data up to 100 mb to be deleted. There is no 
provision for manual adjustment of rejected data; however, opera- 
tional personnel are provided with orbit-by-orbit summaries from 
which they are expected to note trends requiring corrective ac- 
tion. The following tests are performed on the data: : 


Ce Ae I 


a. A superadiabatic temperature lapse rate test is per- 
formed on the tropospheric pressure layers bounded by 
the mandatory and significant pressure levels to a 
insure a stable profile; 


41 


b. A gross error or "neighbor" check is performed on the 
difference between the NMC forecast first guess heights 
and the retrieved heights for all constant pressure 
levels. Values of this difference are calculated for 
all points in the neighborhood (within 500 km). A 
sounding is required to have at least one neighbor or 
it is automatically rejected. The value of the height 
difference for the point is required to agree with the 
average of the height differences for the other points 
in the neighborhood to within: 


+ 200 m with one neighbor, 
+ 100 m with two neighbors, or 
+ 75 m with more than two neighbors. 


c. When retrievals are extended to land areas on an opera- 
tional basis, further location tests will be performed 
to insure that a proper adjustment has been made for the 
terrain elevation in the field of view. 


The above tests are performed on all retrievals; however, the 
initial criteria used in (b) above for the Southern Hemisphere 
may be adjusted to the quality of available first guess profiles 
in this region. 


When a retrieval is rejected, the cause is listed for the infor- 
mation of the operational personnel. At the conclusion of each 
data orbit, the rms temperature difference, E, between the first 
guess and the retrieval is plotted on a latitude-longitude grid, 
where the value of E for a retrieval is given by 


10 A 9.2/2 
E= (1/10) [ = (tT -T 
g=1 


(20) 


an 
where Tg and Tg are the first guess temperatures and retrieval 
temperatures, respectively, for the lowest 10 standard pressure 
levels. 


Histograms of the change made to each first guess temperature 
ere tabulated by mandatory level and by location (the Tropics 
between 18°N and 18°S and Northern Hemisphere and Southern Hem- 
isphere extratropical areas). 
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At the end of each computer processing, coverage charts are 
plotted on polar stereographic maps for the Northern and Southern 
Hemispheres and on a Mercator map for the Tropics. Data rejected 
during the synoptic period are included but flagged in order to 
indicate any consistent geographical bias. 


Additional checks are generated over selected areas or at se- 
lected times to monitor the performance of each VTPR channel, the 
communication links from the data acquisition stations, and the 
individual internal programs used to generate a SIRS sounding. 
Some retrieved soundings are compared with standard radiosonde 
reports, the NMC analyses, and the Scanning Radiometer (SR) data 
as a test for meteorological reasonableness. For example, the 
amount and pressure height of the cloud output are compared with 
the SR data to check the algorithms used in the clear radiance 
program. The 1000 to 300 mb thickness and the contoured height 
fields at 300 mb are also compared with the SR data and the NMC 
forecast and analyses. 


A more complete quality control program is being developed. 
Further details will be published later in a supplement to this 
report. 


5. DATA OUTPUTS AND ARCHIVES 
A. General 


Soundings that pass the quality control tests are archived and 
sent to a number of users. Three tapes are prepared to satisfy 
requirements of various users. One tape is used to send a tele- 
type message; a :ccond tape is sent to the National Meteorological 
Center (NMC); and a third tape is prepared as a data archive. 


B. Teletype Messages 


Teletype messages are formatted to conform to the code used for 
SIRS A and SIRS B (World Meteorological Organization 1972). 
pifferences in the instrument design and the retrieval procedure 
require a change in the interpretation of several of the values 
given. 


A number of adjacent measurements are used to obtain "clear" 
radiances from cloud contaminated values. These "clear" radi- 
ances are used in the retrieval, so cloud conditions do not affect 
the reliability of the data in the way that they affected the 


SIRS soundings. The cloud indicator is always set at zero to in-_ 


dicate a clear sounding, and no cloud information is transmitted. 
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This does not, however, indicate that no clouds are present. 


When the Nimbus 3 and Nimbus 4 satellites were operational, 
retrievals were attempted only in the Northern Hemisphere where 
the NMC forecast was available. The 850-mb forecast height was 
used as a reference level for the height calculations. This pro- 
cedure is currently being used with the VTIPR data for regions 
covered by the NMC forecast. When the 850-mb forecast height is 
not available, the 1000-mb height is set to zero; for these cases, 
the value given as the 300-mb height is actually the 1000- to 
300-mb thickness and the 1000-mb height is zero. 


In summary, north of latitude 21°N all heights are referenced 
to the NMC 850-mb forecast height. South of latitude 18°N, all 
heights are actually thicknesses between the given level and 1000 
mb. Between 18°N and 21°N, either method may be used, and the 
1000 -mb height must be checked to determine which is used for a 
given sounding. 


To make the SIRS code more compatible with the capability of 
the VTPR instrument, the code will be modified in the near future, 
and "clear" radiances will be sent as a Separate transmission. 


C. Output to NMC 


Soundings supplied to NMC are written on two tapes, one for 
each twelve-hour interval (0600-1800 GMT and 1800-0600 GMT). Data 
are output as soon as orbits are processed, and NMC uses the 
tapes as required by their operational schedule. The tapes are 
not available for general use, but the same soundings appear on 
the archival tape. 


D. Archival Tapes 


VTPR data are available in different forms (raw radiances, 
"Clear radiances", and retrieved profiles) at three major points 
in the data processing. These data are written on tape as three 
separate files, each of which consists of a header record and a 
number of data records. The data content and the formats of 
these files are described in appendix III. This tape is trans- 
mitted to the Goddard Institute for Space Studies (GISS) in real 
time. A copy also is sent to the National Climatic Center at 
Asheville, North Carolina. Requests for data and questions con- 
cerning formats should be sent to the National Climatic Center, 
Federal Building, Asheville, North Carolina 28801. 
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APPENDIX I 


Miscellaneous Tables | 


Information contained in tables 2 and 3 applies to the instru- 
ments carried on NOAA 2. Tables 5 and 6 contain information 
about current operational procedures. These tables will be re- 
vised as needed and included in future revisions to this report. 
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Table 1.--VTPR data sequence and spot identifiers, 
(normal sequence) 


Word no. No. of words Position code Description 
1-7 Z 24/25 Primary optics 
8 up 24/25 Not used 
9-16 8 1 Spot 1, filters 1-8 
17-24 8 2 Spot 2, filters 1-8 
25-32 8 3 Spot 3, filters 1-8 
33-40 8 A Spot 4, filters 1-8 
41-48 8 5 Spot 5, filters 1-8 
49-56 8 6 Spot 6, filters 1-8 
57-64 8 7 Spot 7, filters 1-8 
65-72 8 8 Spot 8, filters 1-8 
7 3-80 8 9 Spot 9, filters 1-8 
81-88 8 10 Spot 10, filters 1-8 
89-96 8 1) Spot ll, filters 1-8 
97-104 8 12 Spot 12, filters 1-8 
105-112 8 13 Spot 13, filters 1-8 
113-120 8 14 Spot 14, filters 1-8 
121-128 8 15 Spot 15, filters 1-8 
129-136 8 16 Spot 16, filters 1-8 
- 137-144 8 L7 Spot 17, filters 1-8 
145-152 8 18 Spot 18, filters 1-8 
153-160 8 19 Spot 19, filters 1-8 
161-168 8 20 Spot 20, filters 1-8 
169-176 8 21 Spot 21, filters 1-8 
177-184 8 22 Spot 22, filters 1-8 
185-192 8 23 Spot 23, filters 1-8 
193 1 24/25 Not used 
194-195 2 24/25 Shroud 
196 1 24/25 Not used 
197-199 3 24/25 Secondary optics 
------ eer er ee flyback time - ---+-+-+-+e-+-+-+--- 
200 1 24/25 Primary optics 
1-7 7 24/25 Primary optics 


a 
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Table 2.--Spectral characteristics of VTPR filters, 
1 and 3 


Filter 
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Table 3.--Spectral characteristics of instrument number 1, as 
determined by SEL and manufacturer. 


Central Filter band pass 


wave- Half- Tenth- Hundredth- Maximum 
Filter Laboratory number width width width transmittance 


(cm71) (cm-l) (cm-1) (cm-1) 


1 SEL 667%5 346 8.9 137 232 
BEC 668.0 369 1039 13.4 - 30 
3 SEL 694.8 Load Le Pe, ahead 54 
OCLI 695.1 L265 L749 2162 54 
A SEL 708. 3 10.6 15.0 2138 61 
OCLI 708.2 10.8 1325 20.6 61 
5 SEL I2oa3 11.3 E5.6 20,9 74 
OCLI 725.4 pts PS 16.0 2240 74 
6 SEL 747.4 12.0 16.1 21.4 70 
OCLI 747.4 Lis? 16.8 Ze si3 
5 SEL 532,.80- 15.4 26.4 3143 49 
OCLI 533,65: 15.1 2043 35.8 - 50 
8 SEL 83561 70 15.20 2229 67 
OCLI 835.4 Lok 15:06 2602 - 66 
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Table 4.--VTPR data sequence and spot identifiers (calibration 


321-352 
353-559 
560 
561-567 
568 
569-591 
ays 
593-607 
608 


sequence) 


No. of words Position code 


7 26 Primary optics 
ut 26 zero Offset level 
1 26 Not used 
7 26 lst voltage level 
1 26 Not used 
7 26 2nd voltage level 
i 26 Not used 
fi 26 3rd voltage level 
1 26 Not used 
7 26 4th voltage level 
1 26 Not used 
7 26 5th voltage level 
1 26 Not used 
7 26 6th voltage level 
1 26 Not used 
7 26 7th voltage level 
256 27 Space look 
Filters 1-8 cycling 
32 28 Mirror moving, not used 
207 28 Patch look, filters 12°8 cycling 
ui 28 Shroud 
7 29 Shroud 
1 29 Detector 
23 30 Detector 
i 30 Patch 
as) shal Patch 
1 our Not used 


Description 


Normal scan sequence resumes with spot 1, etc. 


Table 5.--NMC fields used in VTPR retrievals 


Analysis Forecast 
Field interval interval 
Pressure: 
Tropopause 12 hr 12 & 18 hr 
Temperatures: 
Tropopause 12 hr 12 & 18 hr 
1000 to 100 mb 12° hy 12 & 18 hr 
(10 levels) 
70) to.10° mb 24 hr -- 
(4 levels) 
Relative humidity: 
Tropopause i2 -Ae -- 
Boundary layer Le Ar 1i2 & 18 hr 


Heights 
850 mb I2..hx 12 & 18 hr 


Table 6. 
Level Pressure 
(mb ) 

ak -010000 
2 -022509 
| -043472 
4 -075634 
5 - 121989 
6 2185758 
7 ~e1OS7 5 
8 - 379474 
9 - 516882 
10 -686604 
1l -892818 
12 1.43987 1 
13 1.432267 
14 1.774667 
ales) 2el7 Lest 
16 2.628863 
Ly 34150709 
18 3.742052 
19 4.410059 
20 5.158426 
2i 5.993379 
22 6.920666 
23 7.946158 
24 9.075845 
25 10.315836 
26 I1672352 
27 13,151729 
28 14.780413 
29 16.504959 
30 18.392029 
31 20.428391 
32 22.620917 
33 24.976580 
34 27.502455 


Level 


Pressure 


(mb ) 


30°22057.17 
332095637 
36.173585 
39.453026 
427939518 
46.640713 
50.564355 
54.718280 
59.110411 
63.748763 
68.641437 
73.7 90022 
795222593 
84.927709 
90.920415 
97 ..209237 
103.802787 
L102 7097-57 
117.938919 
125.499127 
1333599325 
141.648495 
150.,255758 
159.230270 
168.594.2783 
178. 318103 
188.450141 
198.986865 
2092937822 
Z22leSL2osa 
233.120986 
245.372655 
258.077476 


Level 


--The 100 pressure levels used in VTPR retrievals 


Pressure 


(mb) 


271.245360 
284.886288 
299. 010313 
313.6275596 
328.748216 
344. 382547 
360. 540883 
377 .2383623 
394.471232 
412.264246 
430 .623266 
449 .558960 
469.082061 
489 .203370 
509). 933752 
531 .284137 
953''265520 
575.888960 
599), 1 G>oN9 
623.106566 
647.723168 
673.026699 


699 ..0128533 


725.740107 
753.07 292 
7812338533 
810.248566 
839. 914701 
870. 348681 
901.562308 
933.567446 
966..376016 
1000 .000000 
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APPENDIX II 


Filter Characteristics and CO. Transmittances 


Filter curves have been obtained for the first four VTPR instru- 
ments. These curves were used to calculate CO+ transmittances 
and the resulting weighting functions (derivatives of the trans- 
mittances with respect to a function of pressure). Table 1 
lists the pressures at which numerical values of transmittances 
and weighting functions are given. 


For each instrument (set), twelve pages of data are given in 
the following order: 


filter curves for channels 1 and 2 

filter transmittances for channels 1 and 2 
filter curves for channels 3 and 4 

filter transmittances for channels 3 and 4 
filter curves for channels 5 and 6 

filter transmittances for channels 5 and 6 
filter curves for channels 7 and 8 

filter transmittances for channels 7 and 8 
CO transmittances at 0° zenith angle 

CO, weighting functions at 0° zenith angle 
CO5 transmittances at 23.8° zenith angle 
C09 weighting functions at 23.8° zenith angle. 


Instruments 1 and 3 were used on NOAA 2. 
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Table 1.--The 50 pressure levels for which CO9 transmittances 


Level 


oOo on FF KN 


10 
12 
14 
ue) 
18 
20 
eZ 
24 
26 
28 
30 
32 
34 


and weighting functions were calculated 


Pressure 


(mb) 
.022509 
.075634 
.185758 
. 379474 
. 686604 

1.139871 
1.774667 
2.628863 
3.742652 
5.158426 
6.920666 
9.075845 
112672352 
14.780413 
18. 392029 
22.620917 
27.502455 


Level 


36 
38 
40 
42 
44 
46 
48 
50 
52 
54 
56 
58 
60 
62 
64 
66 
68 


Pressure 


(mb) 
33.093637 
39.453026 
46.640713 
54.718280 
63.748763 
73.796622 
84.927709 
97.209237 

110.709757 
125.499127 
141.648495 
159.230270 
178.318103 
198.986865 
221.312631 
245.372655 
271.245360 


Level 


70 
72 
74 
76 
78 
80 
82 
84 
86 
88 
90 
92 
94 
96 
98 
100 


Pressure 


(mb ) 
299. 010313 
328.748216 
360.540883 
394,.471232 
430 .623266 
469.082061 
509.933752 
553.265520 
599.165579 
€47.723168 
699.028533 
753 17292: 
810.248566 
870.348681 
933.567446 

1000 .000000 
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APPENDIX III 


Archive Formats 


Archival tapes for VTPR data contain earth located raw radian- 
ces, first guess information, and retrievals. Each type of data 
is written as a Separate file. 


The first file contains earth-located raw radiances. These 
radiances have been calibrated - that is converted from counts to 
radiances - but are otherwise unchanged. A description of this 
file is contained in tables 1 and 2. The calibration identi fica- 
tion number referred to in table 1 is now 4. If changes in the 
regression coefficients used for the calibration are required be- 
cause of changes on the spacecraft, this number will change. The 
instrument number referred to in table 2 identifies a scan line 
as: data from the calibration sequence in either instrument (0); 
earth radiances from instrument 1 (1); earth radiances from in- 
stument 2 (2); dual mode earth radiances from the secondary in- 
strument (3). During the rarely-used dual mode, the primary 
instrument controls the data system, and radiances from that in- 
strument are identified normally; radiances from the secondary 
instrument replace those of the primary instrument on alternate 
scans. 

The second file on the tape contains "clear radiances" obtained 
by the CLRAD program, as well as the first guess information that 
was available. A description of this file is contained in tables 
3 and 4, 


The third file on the tape contains retrieved profiles that 
have passed the quality control tests. A description of this 
file is given in tables 5 and 6. Data records on this file dif- 
fer from all other records on the tape in that the data are coded 
and six bits are required to represent each character. Table 7 
shows the octal and binary representations of the characters con- 
tained in the data records. The format used for the data records 
is the same format used for the NMC data, except that twelve 
reports constitute one physical record of 30,600 bits on the NMC 
tape. In order to keep the buffer size small, this record was 
broken down into five physical records of 6,120 bits each for the 
archive tape. 


Table 1.--VTPR archive I format and data description -- 
header record (6480 bits) 


Description 


Satellite number 

Start of data (year) 

Start of data (month) 

Start of data (day) 

Julian day 

Start of data (hour) 

Start of data (minute) 

Start of data (second) 

End of data (hour) 

End of data (minute) 

End of data (second) 

Calibration identification number 

Standard deviations [mw/ (m2sr cm-l) ]x 10 
(channels 1-8 of primary instrument), 
(8 words of fill or channels 1-8 of 
secondary instrument for dual mode 
operation) 

Fill 


ORBITAL PARAMETERS 


Eccentricity x 10° 

Argument of perigee (deg x 103) 
Right ascension (deg x 103) 
Inclination (deg x 103) 
Semi-major axis (km x 103) 
Mean anomaly (deg x 103) 

Epoch year 

Epoch month 

Epoch day 

Epoch hour 

Epoch minute (whole) 

Epoch fraction of minute x 100 
Phi max x 100 

Lambda x 10 

Fill 


Number 
of words 


OR PRP RPRPRPHPHPHPHPHPE 


= 


ee ee 


481 


word 


36 
12 
12 
2 


Bits/ Number 
of bits 


36 
a2 
12 
12 
12 
12 


—— = 
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Table 2,--VTPR archive I format and data description -- 
documentation records (groups of 3 records - 6480 bits per record) 


Description 


2 = 
Radiances [mW/(m~ sr cm 


=) 


] x 20 


(channels 1-8 (spots 1-23 (lines 


1-8) )) 
Latitude (deg + 90) x 10 
East longitude (deg x 10) 
(spot 5 = line 4) 
Latitude (deg + 90) x 10 
East longitude (deg x 10) 
(spot 12 = line 4) 
Latitude (deg + 90) x 10 
East longitude (deg x 10 
(spot 19 = line 4) 
Greenwich Mean Time (GMT) 
Greenwich Mean Time (GMT) 


Greenwich Mean Time (GMT) 


Zenith Angles (deg x 10) 


(hour ) 
(minute) 


(second ) 


Number 
of words 


1472 


(spot 1 = line 1), (spot 1 - line 23), 
(spot 4 = line 5), (spot 4 - line 12), 
(spot 4 = line 19), (spot 8 - line 1), 


(spot 8 - line 23) 
Fill 


Line count 


i. 


1 


Sea surface temperatures [(°K-269.9)x5] 184 


(spots 1-23 (lines 1-8) ) 


her gl 


Instrument number (see text) 


(8 lines) 
Clock 
(spot 1 - line 1) 


Bits/ Number 


word 


2 


12 
i2 


12 
12 


12 
2 


12 


2 


12 


12 


LZ 


12 


28 


36 


of bits 


17664 


i 
2 


12 
i 


12 


L2 


12 


12 


12 


84 


LZ 


1.2 


1472 


24 


36 
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Table 3.--VTPR archive II format and data description -- 
header record (480 bits) 


Number Bits/ Number 


Description of words word of bits 
Start of data (year) id 12 iby 
Start of data (month) ek 12 12 
Start of data (day) ll HZ i2 
Satellite number ul 36 36 
Calibration table number ui 12 12 , 
Instrument number al 12 12 , 
Standard deviations [mw/ (m2 sr cm7!) ]x 100 16 L2 92 
(channels 1-8 of primary instrument) ; : 
(8 zeros or channels 1-8 of secondary 
instrument of dual mode operation) ! 
Start of data (hour) di 12 2 
Start of data (minute) ul 2 2 
Start of data (second) 1 12 iy ia 
End of data (hour) 1 (2 12 : 
End of data (minute) 1 2 2 
End of data (second) 1 2 12 1 ( 


Pat 10 12 120 


Table 4.--VTPR archive II format and data description-- 


documentation records (720 bits) 


Description 


Latitude (deg.+ 90) 


x 10 


West longitude (deg x 10) 


Time of sounding (hour) 


Time of sounding (minute) 


Time of sounding (second) 


ii 
Clear radiances [mW/ (m2 sr cn) 1 x 20 


(channels 1-8) 


Zenith angle (deg x 10) 


Sea surface temperature (°K-269.9) x 5 


First guess values temperatures” (°Cx10) 
(15 standard pressure levels--1000 to 10mb) 


(tropopause temperature) 


Pressure™* (mb x 10) 


Dew point depression! 
(first 10 standard levels) 


Confidence factors 
(channels 1-8) 


First guess 850-mb height”™ (m) 


(x 100) 


obtained from NMC forecast 


Fuld 


*Missing data = 


tMissing data = 


7777 octal 
0000 octal 


Number 
of words 


OF FP FE FP 


1 
a 


16 


10 


word 


TZ 
12 
L2 
2 
12 
12 


2 
12 
12 


TZ 


2 


12 


12 


12 


Bits/ Number 
of bits 


L2 
2 
LZ 
12 
2 
96 


12 
12 
192 


L2 


120 


96 


12 


108 


109 


110 


Table 5.--VTPR archive III format and data description-- 
header record (480 bits) 


Number Bits/ Number 


—— 


Description of words word of bits 
Start of data (year) il 12 We) 
Start of data (month) 1 12 12 

Start of data (day) Ab 12 12 
Satellite number it 36 36 
Calibration table number ut {2 12 
Instrument number L 1Z a2 
Standard pressure levels (mb) 15 2 180 | 
Start of data (hour) a 12 12 | 
Start of data (minute) 1 12 2 

Start of data (second) 1 12 12 | 
End of data (hour) Hi dh2 N2 

End of data (minute) 1 12 12 | 
End of data (second) a 2 2 
Fill 17 12 #32 


hi ig 


Table 6.--VTPR archive III format and data description-- 
documentation records* 


i 


Number Characters/ Bits/ No. of 


Description of words word charac, bits 
cg ee a ee eee arene 
Latitude (deg x 100) i 5 6 30 
West Longitude (deg x 100) i 5 6 30 
Fill (N) 1 1 6 6 


Orbit number for 12-hour 
period (set to 1 at 0600 


and 1800 GMT) 1 1 6 6 
Sounding number 1 2 6 18 
Fil 1 1 6 6 
Time (hours + fraction of 

hours) x 100 I. 4 6 24 
Parameters (data used by 

NMC only) 1 30 6 180 
Data for 15 levels is) 22 6 1980 

Height (m) 

Temperature (°C x 10) 4 


Tg (not given) (°C x 10) 


fad 1 10 
Pil 1 20 6 120 
(End report) 1 10 6 60 
(End cea) variable 10 6 variable 


*A data record consists of 12 soundings or reports and contains 
30600 bits, It is broken up into five physical records of 
6120 bits each, The data on this record are written as an "A" 
format so 6 bits are required for each character (table 7). 


‘when 12 reports are obtained or an orbit is ended, a data 
record is written, The remaining space is filled with "end 
record" so that the five physical records always contain a 
total of 30600 bits, 


TZ 


Table 7.--Six-bit binary code used for data records on the third 
file of the archival tape 


Character Octal representation Binary representation 

0 33 011011 
if 34 011100 
2 35 011101 
| 36 011110 
4 a7, OTT 
5 40 100000 
6 41 100001 
7 42 100010 

8 43 100011 | 

9 44 100100 | 

= 46 100110 | 
=a LOlitT 
Blank 55 LOLOL 
Ze 010010 
E 05 000101 
C 03 000011 
O 17 OOTLIT 
R 22 010010 
D 04 000100 
Pp 20 010000 
a 24 010100 
N 16 001110 
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NESC 


NESC 


NESC 


NESC 


NESC 


NESS 


NESS 


NESS 


NESS 


NESS 


NESS 


NESS 


NESS 


NESS 


51 


ao 


54 


56 


59 


57 


58 


39 


60 


61 


63 


64 
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